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A study is undertaken to develop a methodology for determining the suitability of various 
high-cycle fatigue models for metallic structures subjected to combined thermal-acoustic 
loadings. Two features of this problem differentiate it from the fatigue of structures subject 
to acoustic loading alone. Potentially large mean stresses associated with the thermally pre- 
and post-buckled states require models capable of handling those conditions. Snap-through 
motion between multiple post-buckled equilibrium positions introduces very high 
alternating stress. The thermal-acoustic time history response of a clamped aluminum beam 
structure with geometric and material nonlinearities is determined via numerical simulation. 
A cumulative damage model is employed using a rainflow cycle counting scheme and fatigue 
estimates are made for 2024-T3 aluminum using various non-zero mean fatigue models, 
including Walker, Morrow, Morrow with true fracture strength, and MMPDS. A baseline 
zero-mean model is additionally considered. It is shown that for this material, the Walker 
model produces the most conservative fatigue estimates when the stress response has a 
tensile mean introduced by geometric nonlinearity, but remains in the linear elastic range. 
However, when the loading level is sufficiently high to produce plasticity, the response 
becomes more fully reversed and the baseline, Morrow, and Morrow with true fracture 
strength models produce the most conservative fatigue estimates. 


a 

<r' f , 7 

j ^2? ^3? ^4 


E 

E[P\ 

G 

H 

N. 


Nomenclature 

= coefficient of thermal expansion 

= material fatigue constants (Baseline, Morrow, and Walker Models) 
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= Y oung ’ s modulus 
= expected number of peaks per second 
= shear modulus 

= linear plastic hardening modulus 
= number of cycles to failure 
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P{° min> CT max) 

= joint probability density function (PDF) of stress min and max values 

P{ G ar’ G m) 

= joint probability density function (PDF) of alternating and mean stress values 

G a 

= stress amplitude 

G ar 

= alternating stress 

a eq 

= equivalent stress 

G m 

= mean stress 


= plastic yielding stress 

a TFS 

= true stress at fracture 

R 

= stress ratio (cr min /cr max ) 

RFM ( O’min ’ ^"max ) 

— rainflow matrix 

RFD ( 0"min ? ^"max ) 

= rainflow damage matrix 

P 

= mass density 

T 

= time to failure 

AT 

= temperature increment 


I. Introduction 

F UTURE advanced aircraft and spacecraft structures will be exposed to increasingly severe operating 
environments, including a combination of mechanical, pressure, acoustic, and thermal loads. These loading 
conditions can cause structures to respond in a nonlinear fashion and exhibit complex response characteristics, 
including snap-through behavior, as evidenced by several numerical ’ and experimental " studies. To help ensure 
structural integrity, an evaluation of relevant high cycle fatigue models is required to determine which are most 
suitable over the full range of response conditions. 

Fatigue life prediction for structures subjected solely to high intensity acoustic loading was recently investigated 
by the authors. 6 There it was found that the tensile mean stress associated with the nonlinear bending-membrane 
coupling can significantly reduce the estimated fatigue life. Fatigue life prediction of structures subject to combined 
loadings presents unique challenges because of an even wider range of stress ratios encountered. Consider for 
example, the measured strain response of a thermally buckled plate exposed to acoustic loading. In the absence of 
acoustic loading, the plate assumes one of two thermally buckled equilibrium positions when the temperature 
exceeds the critical value. In this condition, there exists a high static tensile mean strain on one side of the plate and 
a high compressive mean on the other. When low amplitude acoustic loading is applied at this condition, the 
vibratory strain is small compared to the static strain, as shown in the initial portions of Figure 1 and Figure 2. As 
the sound pressure level increases, the magnitude of oscillations increases until an intermittent (Figure 1) or 
persistent (Figure 2) snap-through response occurs. In the intermittent case, either side of the plate experiences both 
tensile and compressive mean strain, whilst the response becomes more fully reversed for the persistent condition. 
As the sound pressure level is subsequently decreased, the response stabilizes about either of the two positions. In 
addition to geometric nonlinearities, material nonlinearity may be present, which changes the global structural 
response and introduces compressive mean stress in locally plastic regions. 

This paper develops a methodology for evaluating the suitability of various high-cycle fatigue models for 
metallic structures subject to thermal-acoustic environments. For this purpose, the thermal-acoustic response and 
fatigue of a flat clamped metallic beam is investigated. This simple structure is considered so that sufficiently long 
stress time histories may be obtained via numerical simulation of the nonlinear dynamic response without being too 
computationally taxing. Both geometric and material nonlinearities are accounted for in the simulation. Four 
distinct fatigue models accounting for the existence of a non-zero mean stress are investigated, namely Walker, 
Morrow, Morrow with true fracture strength, and Metallic Materials Properties Development and Standardization 7 
(MMPDS) models. Additionally, a baseline model neglecting the effect of non-zero mean is considered. The 
fatigue model behaviors as a function of mean stress are examined to understand why a particular model becomes 
more or less conservative than the others under a particular response regime. 
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Figure 1 : Measured strain response of thermally 

2 

buckled plate undergoing intermittent snap-through. 


Figure 2: Measured strain response of thermally 

2 

buckled plate undergoing persistent snap-through. 


II. Thermal-Acoustic Nonlinear Response Analysis 

An aluminum beam structure previously considered 2 served as the basis for the current investigation. The beam 
measured 0.4572 m x 0.0254 m x 2.286 mm (/ x w x h ), and had clamped boundary conditions at both ends, see 
Figure 3. Idealized elastic-plastic material properties with linear strain hardening were used, as shown in Figure 4. 
The elastic-plastic properties were based on empirical data and are given with the elastic modulus, shear modulus, 
plastic modulus, plastic yielding stress, mass density, and coefficient of thermal expansion as: 

E - 73.11GPa, G - 27.59 GPa, // = 417MPa, 

<7 = 400 MPa, p = 2763 a = 22.32 ^4 

P nr C 

It was assumed that the plastic modulus and yield stress are the same in tension and compression. Material 
properties were assumed to be temperature independent because of the modest temperature range considered in this 
paper. Mass proportional damping of magnitude 14.52 s" 1 was used in all the simulations. A failure criterion, e.g. 
ultimate tensile strength, was not imposed on the response modeling. 

The beam response was analyzed with the finite element code ABAQUS. The finite element model consisted of 
144 B21 beam elements 3.175 mm long. The B21 element allows single-plane bending and has one rotational and 
two translational degrees-of-freedom (DoF) at each node. The clamped boundary conditions were modeled by 
constraining the rotational and both translational (transverse and in-plane) DoFs at each end of the beam. 

The ABAQUS/Explicit solution was used with an automatic time step adjustment, referred in ABAQUS as 
‘element-by-elemenf . 8 This approach is known to yield a conservative time step increment. Time histories of the 
displacement, normal stress, and plastic strain, when applicable, were computed for each loading condition. The 
displacements were computed at the mid-span and quarter-span (0.1 143 m) nodes. The stress and the plastic strain 
were computed in proximity to the clamped boundary and the quarter-span location since ABAQUS/Explicit 
computes only element strains and stresses. These locations corresponded to mid-element coordinates of 1.588 mm 
and 0.1159 m, respectively. At each location, the normal stress and plastic strain were recovered at both the upper 
and lower element surfaces. The two locations along the span of the beam were chosen to provide insight into the 
stress response characteristics at differing bending to membrane stress ratios. 

A. Loading 

The beam was simultaneously subjected to a uniformly distributed random pressure loading, p(t ), and a 
uniformly distributed time invariant thermal loading, AT, having zero through-the-thickness variation. Flat, band- 
limited random pressure time histories were generated by summing equal amplitude random phase sine waves at a 
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frequency resolution of 0.61 Hz in the frequency range 0 - 1500 Hz. 9 Three overall sound pressure levels were 
considered: 152, 158, and 170 dB (re: 20 pPa). Each pressure time history was 60 s in duration. 



True Strain s (m/m) 


Figure 4: Elastic-plastic material properties. 

To study the snap-through response, three thermal loading conditions were considered. An ambient temperature 
(AT = 0) case served as an unstressed reference condition. A small temperature increment, AT = 2.78°C (5°F), 
below the critical buckling temperature increment of 3.67°C (6.6°F) was used to investigate the thermally pre- 
buckled response. For both unstressed reference and thermally pre-buckled conditions, a single neutral equilibrium 
position exists. Hence, no snap-through response was expected under any acoustic loading level. 

The thermally post-buckled response was investigated with a temperature increment of 19.44°C (35°F). This 
condition had two non-zero symmetric equilibrium positions, as depicted by the red dashed line in Figure 3. For this 
condition, increasing sound pressure level was previously shown to evolve the response from vibration around one 
of the positions, to intermittent snap-through and finally to persistent snap-through between the two equilibrium 
positions. 2 For each non-zero temperature increment case, the thermal load was applied instantaneously at the 
beginning of each simulation. For the 19.44°C (35°F) temperature increment, application of the thermal loading in 
this manner produced a transient thermal buckling response. 2 

In the remainder of the paper, the combined load case is specified via the shorthand notation x°F/y dB, where x is 
the temperature increment and y is the random pressure level. 

B. Thermal- Acoustic Response Results 

The mid-span transverse displacement response and the total stress response at the clamped end of the thermally 
post-buckled beam (AT = 35°F) are shown in Figure 5 and Figure 6, respectively, to demonstrate the evolution of the 
snap-through response. At the 152 dB level, relatively small amplitude oscillations about one of the thermally 
buckled equilibrium positions is indicated. The initial dynamic thermal buckling transient is visible in the initial 
response. At the 158 dB level, intermittent snap-through response is observed as the mid-span position oscillates 
between the two thermally buckled equilibrium positions. This behavior is akin to the intermittent snap-through 
response experimentally observed in the flat panel. 3 Finally, a persistent snap-through response is shown at the 
170 dB excitation level. A response matrix for the nine loading conditions is shown in Table 1, where N designates 
no snap-through motion, / designates intermittent snap-through motion, and P designates persistent snap-through 
motion. 

As seen in Figure 6, the response at 152 dB and 158 dB was below the plastic yield stress of 400 MPa. As will 
be shown in Section IV, the same was noted in the unstressed reference and pre-buckled conditions at these 
excitation levels. Material nonlinearity was only found to influence the dynamic response when the beam was 
excited at 170 dB, for all three thermal loading conditions. To better understand the effect of plasticity on the 
dynamic response, a simulation was performed at the 35°F/170 dB loading condition using linear elastic material 
properties only. The identical pressure time history was applied as the simulation using elastic-plastic material 
properties. Stress and plastic strain time histories for both simulations are superimposed in Figure 7. The two stress 
responses are identical for the initial 0.15 s or so, deviating thereafter following a stress excursion above the plastic 
yield stress. The point at which they deviate is also marked by a non-zero plastic strain in the elastic-plastic material 
model simulation. Beyond 0.15 s, the elastic material model response continues to oscillate in a symmetric fashion, 
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while the elastic-plastic material model response oscillates mostly about a non-zero mean stress caused by the 
residual plastic strain. The amount of plastic strain varies with time, as subsequent high loading events are 
encountered. The stress response, as shown in Figure 8, exhibits an asymptotic behavior. The maximum stress is 
limited by the amount of peak plastic strain, which increases over time as a result of new high loading events. The 
rate at which these events occur diminishes with time and is governed by the standard deviation of the random 
loading. The onset of plasticity was limited only to a portion of the beam near the clamped ends. While the exact 
extent of the plastic region was not investigated, it was determined that yielding did not occur at the quarter-span. 
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Figure 5: Transverse displacement response at the mid- 
span of the thermally post-buckled beam. 


Figure 6: Total stress response at the clamped end of the 
thermally post-buckled beam. 


Table 1: Response matrix indicating no snap-through 
(N), intermittent snap-through (7), and persistent snap- 
through (P). 
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Figure 7: Stress and plastic strain time history for the 
35°F/170 dB condition at the clamped end. 


The effect of plasticity can also be observed in the stress power spectral density (PSD) at the clamped end. 
Shown in Figure 9 is the stress PSD obtained from the stress response computed using linear elastic and elastic- 
plastic material models. It is seen that the high frequency response characteristics are not significantly altered as a 
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result of plasticity. However, the zero-frequency response shows a substantial increase in compressive mean stress 
due to plasticity. The effect on the fatigue life will be explored further in Section IV. 
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Figure 8: Clamped end stress time history for the 
35°F/170 dB condition. 


Figure 9: Clamped end stress PSD for the 
35°F/170 dB condition. 


III. Random Fatigue Analysis Models 

For a constant amplitude, zero-mean stress response, the stress-life ( S-N) curve relates the cycles to failure, N f , 
to the fully reversed alternating stress, <i ar , via the power law 

a ar = AN b f (1) 


where A = <j' f 2 b , and <j' f and b are material properties determined from tests performed under zero-mean stress. The 
fatigue life can then be determined from 


N f = 


G ar 

A 


( 2 ) 


For conditions having a non-zero mean stress, cr m , a number of models are available, see Dowling. 10 In this 
paper, we consider five models; the baseline zero mean model, Equation (2), two Morrow models, and two Walker 
models. The stress amplitude cr a and mean stress may be expressed as 


o'n =- 


-(i -R) 

qi+tf) 


where the stress ratio R is given by 


For fully reversed loadings, R = - 1, <J ar - <j n 


_ ^min 


( 3 ) 


( 4 ) 


max 

and cr= 0 . 


In the standard Morrow model, the number of cycles to failure is 
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( 5 ) 
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For nonferrous materials, a modification using the true stress at fracture cr TFS , referred to as Morrow TFS in this 
paper, is 


N 


fmTFS 


<J n A 


1 — 


'TFS J 


( 6 ) 


Aluminum 2024-T3 material properties used in the baseline and Morrow models are: 10 cj'j- = 1602 MPa, b = -0.154, 


and <j tfs =610 MPa. 

In the standard Walker model, the number of cycles to failure is 


N f» = 



(7) 


where A w = <j'f w 2 bw and a'f w , b w , and y are material properties determined from tests performed over a range of 

stress ratios R. Material properties used in the Walker model are: 10 <7y W = 1772 MPa, b w = -0.163, and y = 0.460. 

It should be noted that the Walker model material data is not considered valid for R < -2, hence is not calculated 
for stress ratios in that range. 

A modified version of the Walker equation, used in the MMPDS 7 (the replacement for MIL-HDBK-5), is based 
on the stress ratio R and an equivalent stress <j eq , where 


10g(iVy 5 ) = A + A - a 4 ) 

= cr ,„aA ] - R ) h 


( 8 ) 


The term A 4 represents the fatigue limit stress or endurance limit. Section 9. 6. 1.4 of MMPDS 7 gives a thorough 
discussion of Equation (8). It should be noted that the endurance limit is not present in the baseline, Morrow, 
Morrow TFS or standard Walker models. The material constants used in Equation (8), expressed for stress in ksi, 
are: 7 Aj = 11.1, A 2 = -3.97, A 3 = 0.56, and A 4 = 15.8. Because the origin of the MMPDS material property data is 
different than the other models, differences in fatigue estimates may be partially attributable to material properties, 
and not just the model itself. 

The Palmgren-Miner linear cumulative damage rule 11 is typically used for variable amplitude loading and 
assumes that the damage, D, caused by stress cycles for one load amplitude can be calculated and added to damage 
caused by stress cycles in another load amplitude, or 


*>-i 


K 


r(n,\ 


(9) 


where N i are 


the number of cycles, and ( Nf ) 


are the number of cycles to failure for load case i. 


A. Random Fatigue 

For random response, it is convenient to recast Equation (9) in the form 

OO 00 

E[D] = E[P]T | | N p(<J min ,<J mm )d<J min da mm (10) 

— OO — OO 

where E[P ] is the expected number of peaks per second, T is the analysis period, and /?(<T min ,<T ma x ) is the joint 
PDF of stress minimum and maximum values. In terms of amplitude and mean stress 


E[D] = E[P]T J \N f {cr a ,a m r x p(a a ,a m )da a da m . (11) 

— oo 0 
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The dependency of N f on c> m and cr a , or on cr min and cr max , is determined by the particular model used, e.g. 
zero-mean, Morrow, or Walker. 

It was shown 12 that p(c m{n , cr max ) can be estimated from the rainflow matrix (RFM) 

1 OO 00 00 00 

;^XX^ M Lmin>°max)= { J H CT min^max)^min^max ( 12 ) 

RF _oo —00 —00 _Q0 

where N ^ are the number of rainflow cycles. The RFM is a two-dimensional histogram of rainflow ranges. In 

this paper, the RFM was obtained using the Wave Analysis for Fatigue and Oceanography (WAFO) Matlab toolbox 
for analysis of random waves and loads. 13 The archetypal WAFO RFM is shown in Figure 10. Only half of the 
RFM is used because the minimum stress can never be greater than the maximum stress. For example, a stress 
range pair having a zero minimum stress must have a maximum stress greater than zero. Lines are drawn along 
constant values of the stress ratio R to give insight into the distribution of ranges. Constant amplitude conditions, 
when cr max - cr min = constant, are found along diagonals of positive slope. Constant mean stress ranges, when 

°max + ^min = constant, are found along diagonals of negative slope. 



Figure 10: WAFO rain flow matrix. 12 

For a finite record of duration T r , E[P ] « N RF / T r . This is approximate because the discrete RFM disregards 

cycles that have a range less than the bin width. This is considered acceptable given the very low damage from 
these approximately zero amplitude cycles. Substitution of Equation (12) into (10) then yields 


E[D] 



RFM(cJ mm ,a tmx ) 

N f (^min 1 ^"max ) 


J7 00 ' 00 ' 

y X X ’ °max ) 

1 r _oo -oo 


(13) 


where i^FZ>(cr min ,cr max ) is the two-dimensional rainflow damage matrix. Median failure occurs when E[D] = 1 , 
giving the estimated median fatigue life T as 


T = 


tJYL 


RFM((J mm ,cr nmx ) ' 

Nf( a min’^max) 


(14) 


Note that since damage is calculated for each bin of the RFM, the quantities in Equations (13) and (14) vary with 
each cycle. 
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An understanding of how the various models calculate damage for each bin of the RFM is essential for 
interpreting the fatigue results that follow. Plotted in Figure 1 1 - Figure 14 are the cycles to failure for the baseline, 
Morrow, Morrow TFS, and MMPDS models, normalized by the cycles to failure for the standard Walker model. 
The number of cycles to failure were computed directly from Equations (2), (5), (6), (7), and (8), using minimum 
and maximum stress pairs over the range of interest. The figures are therefore not specific to a particular response, 
but are applicable to any condition. The Walker model was selected for the normalization because it was found to be 
the most conservative for the majority of conditions considered. This finding will be discussed in further detail in 
Section IV. Superimposed on each plot is a line for R = -1 . The cutoff in the Walker model below R = -2 is evident 
in each plot. For the baseline model in Figure 1 1, the cycles to failure ratio is largely greater than 1 when the mean 
stress is tensile (to the right of R = -1). A ratio greater than 1 is indicative of a non-conservative estimate relative to 
the Walker model. In the vicinity of R = -1 and for compressive mean stress (to the left of R = -1), the baseline 
model is more conservative than the Walker model. The ratio for the Morrow model is shown in Figure 12. Like 
the baseline model, the Morrow model is non-conservative for a tensile mean stress and more conservative for a 
compressive mean stress. The Morrow model is seen to be more conservative than the baseline model in the 
presence of a tensile mean stress. The Morrow TFS ratio is shown in Figure 13 to exhibit similar characteristics as 
the Morrow model in the highlighted region. This region encompasses most of the cases considered, except for the 
170 dB conditions at the clamped end. For the high stress ranges associated with the latter conditions, the Morrow 
TFS model becomes more conservative than the Walker model. Particularly noteworthy is that the Morrow TFS 
model is more conservative than the Morrow and baseline models for high stress amplitudes along R = -1. Finally, 
the MMPDS ratio is shown in Figure 14. For sufficiently high stress amplitudes, the MMPDS model becomes much 
less conservative than the Walker model along R = -1. The peculiar behavior for small stress amplitudes is 
attributable to the endurance limit present in the MMPDS model. 
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Figure 11: Baseline model cycles to failure normalized by Walker model cycles to failure (Nj- / N ^ ). 
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Figure 12: Morrow model cycles to failure normalized by Walker model cycles to failure ( N^ m / N ^ w ). 
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Figure 13: Morrow TFS model cycles to failure normalized by Walker model cycles to failure ( N^ mTFS / N ). 
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Figure 14: MMPDS model cycles to failure normalized by Walker model cycles to failure ( 7Vy 5 / N ^ ). 


IV. Fatigue Results 

Studies were first conducted to verify that a sufficient quantity of data were available to generate converged 
fatigue estimates. Fatigue estimates were then computed using the aforementioned models to determine which 
tended to be most and least conservative under each response condition, i.e., unstressed reference, pre-buckled, and 
thermally post-buckled. Finally, the effect of localized material nonlinearity was studied to quantify which factor 
affects fatigue life more; the plastic zone itself or changes in the structural dynamics due to a portion of the structure 
having undergone plastic deformation. 


A. Convergence Study 

The goal of the analysis was to estimate time (or number of cycles) to failure with an allowable random error, 
err N j - , of approximately 10%. The power law relationship between stress and number of cycles to failure, Equation 

(1), was used to estimate the random error in the standard deviation of stress response err a 

1 + err G = (l + err Nf ) (15) 

which for this study is approximately 2%. This value of err G was used to set the simulation duration to 60 s. 14 

The initial 1 s segment containing the startup transient was removed from each simulation block, leaving a 59 s 
block of fully developed response for use in ensuing fatigue analyses. The convergence was studied by dividing the 
simulation results into sub-blocks and computing the time to failure estimates as each additional sub-block was 
included in the ensemble. In seven of the nine conditions studied, a single 59 s block was found to be sufficient for 
converged fatigue estimates. A typical plot of this convergence study is shown in Figure 15. Note that the order of 
the time to failure estimates, based on the different damage models, does not change and that the desired -10% error 
for time to failure has been achieved. For the case shown, the stress was below the endurance limit of the MMPDS 
model, so those results are not shown. 

The two cases, 35°F/152 dB and 35°F/158 dB, did not tend to converge for the baseline simulation period of 
59 s, so nine additional simulations were run to determine the time required to reach a converged estimate of time to 
failure. A convergence plot for the 35°F/158 dB case is shown in Figure 16. In both the 35°F/152 dB and 
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35°F/158 dB cases it took approximately 200 s to converge to a stable estimate with the desired -10% error for time 
to failure. Again, the order of estimate based on damage model did not change as the time to failure converges. 




Simulation Time (s) Simulation Time (s) 


Figure 15: Convergence study for a single 59 s Figure 16: Convergence study for ten 59 s simulations 

simulation at the 5°F/152 dB loading condition. at the 35°F/158 dB loading condition. 


B. Unstressed Reference and Pre-Buckled Conditions 

Fatigue estimates at the clamped end are first considered, as the stresses there are highest and will dictate the 
fatigue life of the structure as a whole. The 0°F/152dB loading condition exhibits the effect of geometric 
nonlinearity without the effect of material nonlinearity. The RFM is shown in Figure 17. For small stress 
amplitudes, the majority of cycle counts are approximately evenly distributed about R = - 1. However higher stress 
amplitudes exhibit a tensile mean which is attributable to membrane stretching as the structure undergoes large 
deformation. The effect of this tensile mean on the accumulated damage is apparent in Figure 18 and Figure 19, for 
the Walker and Morrow TFS models, respectively. As shown in Figure 11 - Figure 13, the Walker model 
accumulates more damage for high stress cycles with tensile mean than the Morrow TFS, Morrow, and baseline 
models in this stress range. The stress range for this condition is below the endurance limit of the MMPDS model. 
The fact that damage accumulation for the Walker model is cutoff for R < -2 is of little consequence, as the damage 
associated with these cycles is small. The fatigue life estimates for this condition are shown in Figure 20. The 
Walker model is clearly the most conservative. However, it should be noted that all estimates are within a factor of 
three. In this and similar figures that follow, dashed lines were used to highlight the trends. Interpolation along 
these lines to establish intermediate values is not valid. 



Figure 17: RFM for the 0°F/152 dB condition. 


Figure 18: Walker model RFD for 0°F/152 dB. 
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Figure 19: Morrow TFS model RFD for 0°F/152 dB. 
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Figure 20: Clamped end fatigue estimates at 152 dB. 


As the acoustic excitation level is increased from 152 dB to 158 dB, the effect of geometric nonlinearity is 
increased, while the stress remains below the plastic limit. The corresponding RFM is shown in Figure 21. For the 
reasons cited for the 0°F/152 dB conditions, the Walker model remains the most conservative of all models, as 
shown in Figure 22. The stresses for the 0°F/158dB condition exceed the endurance limit of the MMPDS model. 
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Figure 21: RFM for the 0°F/158dB condition. 
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Figure 22: Clamped end fatigue estimates at 158dB. 


As the acoustic excitation level is further increased to 170 dB, stresses are sufficiently high to cause the material 
to yield in the clamped region. As previously noted in Figure 9, plasticity results in a compressive mean stress, 
which, when combined with the tensile membrane component, shifts the RFM toward the R = -1 main diagonal, as 
shown in Figure 23. While the low stress cycles now exhibit a compressive mean, the high stress cycles lie close to 
R = - 1, due to the combined effect of material and geometric nonlinearity. The majority of accumulated damage 
consequently lies close or on R = -1 . A comparison of the Walker and Morrow TFS RFD matrices, in Figure 24 and 
Figure 25, respectively, shows that the Morrow TFS model accumulate more damage than the Walker model. 
Recalling Figure 11 and Figure 12, this is also true of the baseline and Morrow. Only this MMPDS model is less 
conservative in this stress region. The corresponding fatigue life estimates are shown in Figure 26. Here the fatigue 
life estimates are within a factor of two, reflecting the fact that the difference in damage accumulation is less 
prevalent near R = -1 . 
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Figure 23: RFM for the 0°F/170 dB condition. 




Figure 25: Morrow TFS model RFD for 0°F/170 dB. 




Figure 24: Walker model RFD for 0°F/170 dB. 
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Figure 26: Clamped end fatigue estimates at 170 dB. 


Finally, in the pre-buckled configuration, there exists a compressive mean stress of approximately 57 MPa due to 
the thermal loading alone. Due to the dynamic loading, the structure responds in a geometrically nonlinear fashion 
and the mean stress diminishes. Like the plastic condition, this shifts the RFM to the left. However, the amplitude 
is not high enough to cause the peak stress region to move sufficiently close to R = -1 to make a significant 
difference in damage accumulation. Hence, the fatigue life estimates are similar to those noted in the unstressed 
reference condition, with Walker the most conservative in the linear elastic range, and Morrow, Morrow TFS, and 
baseline models more conservative in the plastic range. The fatigue life is reduced for this condition relative to the 
unstressed reference condition as the compressive stress causes the beam to become less stiff, resulting in a higher 
stress response for a given acoustic excitation level, see Figure 20, Figure 22, and Figure 26. 


C. Thermally Post-Buckled Condition 

The thermally post-buckled condition presents more interesting dynamic response behavior, resulting in either 
large amplitude oscillations about one of the equilibrium conditions, intermittent snap-through between equilibrium 
positions, or persistent snap-through between positions, as previously shown in Figure 5 and Figure 6. 

Consider first the 35°F/152 dB condition which results in large amplitude oscillations about one of the 
equilibrium positions. This is the only condition studied in which the top and bottom surfaces of the beam have 
differing time averaged stress responses. The RFM for the surface exhibiting a large tensile mean stress is shown in 
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Figure 27. The opposing surface is in compression and will not dictate the fatigue life. Hence, it is not considered 
in the following. The mean stress puts the highest stress cycles in a region in which the Walker model is most 
conservative. This is clearly evident when comparing the Walker model RFD in Figure 28 with the Morrow TFS 
model RFD in Figure 29. The fatigue life estimates are shown in Figure 20 and indicate that the Walker model 
produces the most conservative fatigue estimates. 



-200 -150 -100 -50 0 50 100 150 200 
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Figure 27: RFM for the 35°F/152 dB condition. 


Figure 28: Walker model RFD for 35°F/152 dB. 


As the sound pressure level increases to 158 dB, the response enters into an intermittent snap-through motion 
between two symmetric equilibrium positions. This produces an RFM resembling a wishbone, as shown in Figure 
30. Note that a single 59 s ensemble was used to generate this figure and the corresponding RFDs for clarity, but 
fatigue life estimates were made with 10 ensembles, as indicated in Section IV. A. 
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Figure 29: Morrow TFS model RFD for 35°F/152 dB. Figure 30: RFM for the 35°F/158 dB condition. 


Damage accumulation for the Walker and Morrow TFS models are shown in Figure 31 and Figure 32, 
respectively. Common to both models is a lack of damage accumulation associated with the lower equilibrium 
position (for the upper surface of the beam). Here the upper surface is in compression. The Walker model is seen to 
accumulate more damage for cycles oscillating about the upper equilibrium position than the Morrow TFS model, as 
was found to be the case for the 35°F/152 dB condition. The most damaging events are associated with the large 
stress ranges as the structure snaps-through from one equilibrium position to the other. However, as this occurs only 
intermittently, most of the accumulated damage is still associated with oscillations about the upper position for the 
Walker model. Because of this, the Walker fatigue life estimate is still the most conservative of all the models, as 
shown in Figure 22. 
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Figure 31 : Walker model RFD for 35°F/158dB Figure 32: Morrow TFS model RFD at 35°F/158dB. 


Lastly, consider the 35°F/170 dB condition. Here, the stress ranges which lie between the equilibrium positions 
are filled in by persistent snap-through motion, as shown in Figure 33. The RFM is much broader than either the 
0°F/170 dB or 35°F/170 dB conditions because of the oscillations about the two equilibrium conditions. However, 
because of plasticity, the maximum stress is only slightly higher than either of those conditions. The broadness of 
the RFM is not significant however, as the majority of damage is along the fully reversed bending line R = - 1. Here, 
as in previous 170 dB conditions, the damage accumulated by the Morrow TFS model is greater than that of the 
Walker model, as shown in Figure 34 and Figure 35. The fatigue estimates for all models are shown in Figure 26. 



Min Stress (MPa) 



Figure 33: RFM for the 35°F/170 dB condition. 



Figure 34: Walker model RFD for 35°F/170 dB. 


Fatigue estimates made for the quarter-span location at shown in Figure 36, Figure 37, and Figure 38 for 152 dB, 
158 dB, and 170 dB loadings, respectively. Because the stress is lower at the quarter-span than at the clamped end, 
all loading conditions indicate longer fatigue lives than the clamped end counterparts. In all but the 170 dB case, the 
stress is below the MMPDS endurance limit, and therefore MMPDS estimates are not available for the majority of 
conditions. Because of geometric nonlinearity, the tensile mean makes the Walker model produce the most 
conservative estimates. Unlike the clamped end, the quarter-span does not yield under any of the 170 dB conditions. 
Therefore, there is no compressive mean from material nonlinearity to offset the tensile mean from the geometric 
nonlinearity. A slight rise is seen in the Walker model estimate at the 5°F/152 dB condition. For this condition, 
there is a compressive mean stress. As the compressive stress increases, all models except the baseline model 
accumulate less damage for a given amplitude. The rate at which this trend occurs is much higher for the Walker 
model than for the other models. It is seen that the differences in fatigue life between the various models at the 
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quarter- span are smaller than the differences observed at the clamped end for all three conditions with 35°F 
temperature increments. 



Figure 35: Morrow TFS model RFD for 35°F/170 dB. 
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Figure 36: Quarter-span fatigue estimates at 152 dB. 
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Figure 37: Quarter-span fatigue estimates at 158 dB. Figure 38: Quarter-span fatigue estimates at 170 dB. 


D. Effect of Localized Material Nonlinearity 

The ratios between fatigue estimates obtained with the elastic material model and estimates obtained with the 
elastic-plastic material model are presented in Table 2 for all three 170 dB conditions at the quarter- span and 
clamped end. It is seen that the elastic-plastic material model generally yields longer fatigue estimates than the fully 
elastic model. Only two estimates made with the baseline model (5°F and 35°F) indicate a slightly longer life for 
the elastic material model. Thus, material nonlinearity has a beneficial effect of limiting high amplitude stress 
cycles, albeit with the detrimental effect of permanent plastic deformation. Because the quarter-span location 
remained in the elastic response regime, material nonlinearity did not directly limit the stress in this location. 
Instead, the benefit was a result of a reduction in the dynamic response due to yielding nearer to the clamped end. 
This conclusion is borne out in Table 2, which shows a greater reduction in elastic model fatigue life at the clamped 
end, where both local plasticity and global dynamics are contributors, than at the quarter-span where only global 
dynamics contribute. Compare, for example, the Walker model at ambient conditions which indicates at clamped 
end fatigue life ratio of 0.360, while the quarter-span ratio is 0.799. 
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Table 2: Fatigue estimates obtained with elastic material model at 170 dB. 


AT and 
Location 

Fatigue Life Ratio 

Elastic Material Model / Elastic-Plastic Material Model 



Baseline 

Walker 

Morrow 

Morrow TFS 

MMPDS 

0°F 

CL 

0.638 

0.360 

0.436 

0.219 

0.519 

Q 

0.815 

0.799 

0.789 

0.741 

0.562 

5°F 

CL 

0.935 

0.557 

0.688 

0.411 

0.684 

Q 

1.013 

0.997 

0.991 

0.952 

0.871 

35°F 

CL 

0.608 

0.410 

0.450 

0.240 

0.396 

Q 

1.025 

0.903 

0.992 

0.935 

0.877 


V. Conclusions 

Five distinct high-cycle fatigue models were considered in the analysis of an aluminum beam subjected to 
combined thermal-acoustic loadings. Four of the models accounted for the presence of mean stress, and one 
neglected this effect. The structural dynamic response regimes included an unstressed reference condition, a pre- 
buckled condition, and a thermally post-buckled condition which permitted snap-through behavior. A broadly 
applicable methodology was developed to determine the most conservative fatigue model through an understanding 
of how each model accumulates damage as a function of the mean stress. For the specific problem considered, it 
was found that the Walker model provided the most conservative fatigue estimate when the stress response had a 
tensile mean due to the effect of geometric nonlinearity alone. For sufficiently high loading conditions, the Morrow 
and Morrow TFS models were found to produce the most conservative estimates in locally plastic regions. These 
conditions were characterized by the most damaging stress cycles being fully reversed. Although not generally 
recommended for nonlinear response problems, the zero-mean baseline model was coincidentally found to produce 
conservative estimates for the high loading conditions. Local plasticity was determined to be beneficial to fatigue 
life relative to the elastic material model. Additional work is needed to generalize the specific conclusions for other 
material systems, structural configurations, and loading conditions. 
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